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Abstract 

A mixture of spin-polarized light and heavy fermionic atoms on a finite size 2D optical lattice is 
considered at various temperatures and values of the coupling between the two atomic species. In 
the case, where the heavy atoms are immobile in comparison to the light atoms, this system can be 
seen as a correlated binary alloy related to the Falicov-Kimball model. The heavy atoms represent 
a scattering environment for the light atoms. The distributions of the binary alloy are discussed 
in terms of strong- and weak-coupling expansions. We further present numerical results for the 
intermediate interaction regime and for the density of states of the light particles. The numerical 
approach is based on a combination of a Monte-Carlo simulation and an exact diagonalization 
method. We find that the scattering by the correlated heavy atoms can open a gap in the spectrum 
of the light atoms, either for strong interaction or small temperatures. 
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INTRODUCTION 



Recent experimental progress in preparing and measuring clouds of ultracold atoms in 
magnetic traps has opened a new way of studying bosonic and fermionic many-particle 
quantum states. Among them are condensed and Mott-insulating states of bosons in optical 
lattices (see e. g. Q, ^, Q]). In comparison with similar studies in solid-state physics, 
atomic clouds enable us to design new many-particle systems by mixing different types of 
atoms P], 0, 0. These mixtures can form new quantum states due to the competition 
between the different types of atoms. In this paper we propose a mixture of light fermionic 
(e.g. ^Li) atoms and heavy fermionic (e.g. ^°K) atoms, and study its low-temperature 
behavior in a (finite) optical lattice. We assume that the cloud of this mixture is prepared 
in a magnetic trap such that the atoms are spin polarized. 

The difference of the masses of the two types of atoms implies two different and well- 
separated time scales for their tunneling processes through the optical lattice. The relatively 
fast tunneling processes of the light atoms sets the relevant scale for the dynamics of the 
mixture. In contrast, the relatively slow tunneling processes of the heavy atoms are dy- 
namically irrelevant and lead only to statistical fluctuations which drive the system towards 
equilibrium. The latter will be discussed by the fact that the heavy particles form Ising-like 
(para-, ferro- and antiferromagnetic) states. They provide a scattering environment for the 
light atoms. Formally, this physical picture leads to the Falicov-Kimball model which has 
been used to describe complex solid-state systems say , lUj, A numerical study of 



the two-dimensional Falicov-Kimball model, based on exact diagonalization, has revealed the 
possibility of a discontineous transition between ordered and disordered phases In the 
following we will use analytic methods as well as a combination of exact diagonalization and 
Monte-Carlo simulations to study large two-dimensional clusters for a better understanding 
of the underlying physics. 

The paper is organized as follows: In Sect. 2 the model is briefly discussed, based on a 
functional- integral representation. A mapping to a binary-alloy model is described in Sect. 
3. This gives the foundation for our analytic treatment, based on weak- and strong-coupling 
expansions, and for the construction of our numerical method. The numerical approach is 
used to evaluate the distribution of the heavy atoms and the density of states of the light 
atoms. 
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THE MODEL 



The atomic degrees of freedom are given in second quantization by local creation and 
annihilation operators. In the case of spin-polarized fermionic atoms we use cj (c^) and 
(fr) as the creation (annihilation) operators for the light and the heavy fermionic atoms, 
respectively, where r denotes the coordinates of the site in the optical lattice. The light 
atoms can tunnel with tunneling rate i, and we assume that the tunneling rate of the heavy 
atoms is so small that we can neglect it. Moreover, there is only a local interaction between 
the atoms in the optical lattice, i.e. only atoms in the same potential well notice each other. 
Since the atoms are spin-polarized fermions, there can be at most one atom per sort in each 
potential well, thanks to Pauli's principle. The interaction strength between light and heavy 
atoms is U. This allows us to write the many-particle Hamiltonian as 



where (r, r') means pairs of nearest-neighbor lattice sites. We have assumed the same chem- 
ical potential n for both types of atoms. This may not be very general but will serve for the 
purpose of studying competing quantum phases in the atomic mixture. The model defined 
in Eq. is also known as the spinless Falicov-Kimball model 

mum 

. It is known 

to describe ordered phases and phase transitions for correlated electronic systems and was 
recently investigated intensively in the limit of infinite dimensions j^. We will study this 
model in the following for a finite lattice, using a correlated binary-alloy (CBA) approach. 

Functional-Integral Representation 

A grand-canonical ensemble of a mixture of light and heavy fermionic atoms at the inverse 
temperature (3 = l/ksT can be defined by the partition function 



where Tr is the trace with respect to all many-particle states in the optical lattice. The 
Green's function for the propagation of a light particle in the background formed by the 
heavy atoms in imaginary time t is 




(1) 



(r,r'> r 



Z = Tre-^^, 



G(r,t;r',0) = ^Tr [e-^^-^^^ce-^^cJ," . 
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These expressions can also be written in terms of a functional integral on a Grassmann 
algebra 13|. For the latter the integration over a Grassmann field '^ai^, t) and its conjugate 
'^^{r,t) (cr = c, /) is given as a linear mapping from a Grassmann algebra to the complex 
numbers. At a space-time point (r, t) we have for integers k,l > 



[^^{r,t)f[^^{r,t)yd^^{r,t)d^^{r,t) = 4,i<^U- 
The partition function Z of the grand-canonical ensemble then reads 

Z = J exp{-S)V[^j,^,] (2) 

with the action 

S=Y^ ^^(r, t) [*.(r, t) - ^,{r, t - A)] + A ^ H[^„{r, t), ^^(r, t - A)] (3) 

r,t,a t 

and the product measure 

r,t,a- 

The Green's function of light atoms is 

The discrete time is used with t = A,2A, (3, implying that the limit A ^ has to 
be taken in the end and (3' = (3 /A is the number of time steps. '^cr{i^,t) and \E'o-('", i^) 
are independent Grassmann fields which satisfy antiperiodic boundary conditions in time 
\l/o-(r, /? + A) = — \I'o-(r, A) and ^^{r,P + A) = — ^cr(r, A). For the subsequent calculations 
it is convenient to rename \I^o-(r, t) '^„{r,t + A) because then the Grassmann field appears 
with the same time in the Hamiltonian of the action Q. 



THE CORRELATED BINARY ALLOY 



The functional integration in Z and G can be performed in several steps, beginning 
with the integration of the heavy atomic field introducing the Ising spins and finally 
integrating the light atomic field \l/c [15]. The details of this procedure are described in 
Appendix A. As a result we obtain for the partition function 

Z= Z{{Sr}) 
{S{r)} 
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with 

Z{{Sr}) = /i^'S.(i+5M)/2det(-5i + Ji + i- (f/72/i)(l + S)). (5) 

The parameters are U' = AU, ft = 1 + A/i, and i is the tunnehng term muhiphed by 
A. Moreover, dt is the time-shift operator. The Ising spin S{r) corresponds with a local 
occupation number nf{r) of the heavy atoms as 

nfir) = [l + Sir)]/2. 

The Green's function is now an averaged resolvent 

G = {{-dt + fi + i- {U'/2fi){l + 5))-^)ising, (6) 

where the average (...)ising is taken with respect to the distribution 

p[{S{r)}) = ^^^^'-^^ (7) 
Y.{S{r)} Zi{Sr}) 

The partition function can also be written as 

{Sir)} 

The distribution is not Z2 invariant (i.e. invariant under a change S{r) —S{r)), except 
for a half-filled lattice (i.e. fi = U/2). 

The representation of the Green's function in Eq. ® is our main analytic result. It means 
that the light particles tunnel through the optical lattice where they are scattered by the 
heavy particles. The distribution of the heavy particles is given by the distribution shown 
in Eq. (jZj). The latter depends on the temperature but also on the parameters of the 
light particles like t and the coupling U between the light and the heavy particles. This 
reflects the intimate relationship between the two types of particles. In other words, the 
light particles move in a random potential formed by the heavy particles. This randomness, 
formally expressed by the Ising spins, is correlated and can be called correlated binary alloy 
(CBA). There is a correlation length which diverges at the phase transitions of the Ising 
system. In the subsequent investigation we will study these phases and their implications 
for the properties of the light particles. 

The symmetric matrix fi—t+{U'/2fi){l + S) can be diagonalized with eigenvalues 1 — AAj. 
Then the determinant in Eq. (jH)) is for A ~ 

det[l + {fx + i- {U'/2ji){l + S)Y'] ~ + ^"''^O- 

i 
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Since the matrix depends on the fluctuating Ising spins, it is diflicult to determine the eigen- 
values. One way to get an idea about the physics of this model is to study the asymptotic 
regimes of strong and weak coupling, another one is to use a numerical diagonalization 
procedure. Both approaches shall be applied subsequently. 



Approximations of the CBA Distribution 

The distribution was studied in the case of strong coupling (i.e. the tunneling (or U~^) 



expansion) in a number of papers [Ij, ll5[. It leads at half-fllling to an Ising model with 
Z2 symmetry. In the following we study the CBA distribution in weak and strong-coupling 
approximations as well as numerically by a Monte-Carlo simulation. The density of states 
of the light atoms are evaluated by a numerical procedure. 



System without Tunneling 

The absence of the Z2-symmetry can be observed already in the absence of tunneling. 
Then we have in the limit A — 

^P^l{l+S{r))/2 _^ g/3[;.+ (;.-C/)(l+5(r))/2] 

PoiiSir)}) = n 1 + 2e/^M + e^i^.-u) 

which is Z2 invariant only for /i = U/2. The average spin is shown in Fig. ^ and its 
asymptotic low-temperature behavior is 

-1 < 



0< fi<U 
1/3 fx = U 

1 U < fx 



Thus only for < fi < U the Ising state is paramagnetic. The other regimes are ferro- 
magnetic. In terms of the conflgurations of the heavy atoms there are no heavy atoms for 
/i < and a heavy atom at each site for fi > U. In the intermediate regime < /i < f/ we 
anticipate that the coupling of neighboring Ising spins, caused by a non-zero tunneling rate 
i, will lead to ordered Ising spins at low temperatures. 
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FIG. 1: Average spin for a system without tunneling: U = 1 and /? = 10. 
Tunneling Expansion 

The effect of a weak tunnehng rate can be evaluated in terms of a perturbation theory 
with respect to t. Moreover, we take the limit A ^ and consider the asymptotic regime 
of low temperatures (i.e. /3 ~ oo). If we include tunneling terms up to order we get 

for ij, = U/2 the Ising model with nearest-neighbor spin interaction: 

Pi^Eir,r')Sir)S{r')+o{P) 



Ps{{S{r)}) 



exp 



(9) 



E{5W=±i} exp [ - P§j Eiry) S{r)S{r') + o{t^) 
This model has an antiferromagnetic low-temperature phase. 

The spin-spin coupling is exponentially small in (3 for fi < and fi > U but of order P/U 
for < < [/. In particular, we can distinguish three different regimes: 

< : Ps{{S}) « exp [ - ^ S{r) + o(P) 



<//<[/ : Ps{{S}) oc exp 
and 



-13U/2 



smh[/3{^i-U/2)]Y,S{r)-/3^ ^ S{r)S{r')+o{^) 
l^-U 



<r,r'> 



U<i,: Ps{{S}) oc exp (3^^^S{r) + o{P' 

^ r 

Besides the two ferromagnetic regimes for // < and /i > U we have the intermediate regime 
< II < U with antiferromagnetic ordering. There is an exponentially small magnetic field 



for /i 7^ f//2 which breaks the symmetry. As we approach fi = or fi = U the magnetic 
field becomes larger. There is a first order transition from the antiferro- to a ferromagnetic 
phase when the magnetic field starts to dominate the spin-spin interaction. This can be 
seen in a simple mean-field approximation. 



The Weak- coupling Limit 

In the case of weak coupling we can perform an expansion in terms of the coupling 
parameter U. This gives in leading order an uncorrelated binary alloy: 

^f3i^l-Ug)S{^)/2 
PUiSir)}) = n y Pii,-Ug)Sir)/2 

where e{k) is the dispersion of the tunneling term and 

^7(/i)=/e(e(A:)+/x)^, 
i.e. < (? < 1. Thus the Ising groundstate for weak coupling is ferromagnetic. In particular, 

(5)=tanh(/3(/i-[/(7(/i))/2). 



The Density of States 

The density of states (DOS) for the light particles can be obtained from the diagonal 
elements of the Green's function in Eq. ©. Its qualitative behavior depends strongly on 
the state of the heavy particles: In the case of a ferromagnetic state the DOS shows a single 
band, for the antiferromagnetic state it has a gap. For a paramagnetic state of the heavy 
particles the form of the DOS is less obvious. We have calculated the DOS numerically 
for a 18 X 18 square lattice with open boundaries. For this purpose we have generated 
configurations of the Ising spins according to the distribution function in Eq. ((Tj), using the 
Metropolis algorithm. Typical configurations with large statistical weight are shown in Figs. 
El -El for different values of the physical parameters [/, /i, and (3. 

For a given configuration {5'(r)} of Ising spins the Hamiltonian 



h = i- 
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FIG. 2: Paramagnetic Ising-spin configuration for t = 1, U = 3 , = U/2, and (3 = 3. White 
(black) squares refer to sites (un)occupied with a heavy atom. 




FIG. 3: Mixture of para- and antiferromagnetic Ising-spin textures for t = 1, U = 3 , n = U/2, 
and (3 = 7 

is diagonalized. From the eigenvalues Xk{{S{r)}) the DOS of h is calculated as 

1 ^ 

D{E,{S{r)}) = —j:HE-Xk{{S{r)})) , (10) 

k=l 

where is the number of lattice sites. Finally, the DOS related to the Green's function in 
Eq. (jH)) is determined by averaging over L = 100 spin configurations: 

D{E) = j Y: D{E,{S{r)}) . (11) 

{S{r)} 



9 



FIG. 4: Antiferromagnetic Ising-spin configuration for t = 1, U = 3 , fi = U /2, and /? = 14 



I 




FIG. 5: Ferromagnetic Ising-spin configuration for t = 1, U = 8 , ^ = 0.8U, and /3 = 14 
In the following the hopping rate is set to t = 1. 

Fig. ini shows the DOS of the light particles for U = 3 and half filling (i. e. /i = U/2) at 
different temperatures. For small /3, i. e. high temperatures, the system shows a gapless 
metallic band, which is symmetric around the Fermi level, and the Ising spins form a para- 
magnetic state. The DOS is slightly suppressed at the band center due to the interaction 
between the light and the heavy atoms. When the temperature is decreased, the Ising spins 
start to order antiferromagnetically, i. e. the heavy atoms create a chessboard-like phase 
with empty sites. This is accompanied by the formation of a gap around the Fermi level 
and a strong enhancement of the DOS at the inner band edges. Very similar results were 
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FIG. 6: DOS for U = 3 , /i = U/2 (half filling) and different temperatures. First row: /3 = 3 , /3 = 
7 . Second row: (3=10, (3 = U. 

found in a dynamical cluster approximation on an 8 x 8 cluster 

The high-temperature regime of the system at half filling is depicted in Fig. [7|for various 
interaction strengths U . For small interaction the DOS shows a metallic band and is peaked 
around the Fermi level. For increasing U this peak gets suppressed and a band splitting to 
two symmetric bands occurs. The spectral weight within these subbands is highest at their 
center. A further increase of U leads to a shift of the lower and upper band to lower and 
higher energies, respectively. 

Figure |H1 shows the DOS in the low-temperature regime for two values of the interaction 
strength {U = 3, solid and U = 8, dashed) and different values of the chemical potential. 
For the latter the distribution ((7j) has no Z2 symmetry, i.e. is not invariant under a global 
spin flip S — > —S. Near half filling the heavy atoms order in a chessboard configuration. 
This behavior is stablized for larger deviations fom /i = f//2 when the interaction strength 
is increased. As /i deviates even further from U/2 the Ising spins start to order ferromag- 
netically. The spectral weight locally shifts to the center of each subband and globally shifts 
from the lower to the upper band. For the completely ordered Ising spins the lower band 
disappears and the system has only one band. 
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FIG. 7: DOS for /3 = 3 at various values of the interaction U and half-filling = U/2). First row: 
U = l, U = 4,. Second row: U = 6, U = 8. 

CONCLUSIONS 

A mixture of light and heavy fermionic atoms is studied as a system in which the light 
atoms live in a correlated disordered environment. This environment is formed by the 
heavy atoms. The disorder is given by fluctuating Ising spins with a complex temperature- 
dependent distribution. This distribution is Z2 (spin-flip) invariant only at half filling (i.e. 
H = U/2) but has a broken spin-flip symmetry for n 7^ U/2. This symmetry breaking 
favors an Ising spin Sr — —1 at low density (i.e. // < 0) and Sr — 1 Bi high density (i.e. 
Ii> U). There is an intermediate regime where an antiferromagentic (staggered) Ising-spin 
configuration is favored. Scattering on these configurations opens a gap in the band of the 
hght atoms. 
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A. Appendix 

A.l. Integration of the Heavy Atoms 

It is possible to integrate out the field of the heavy atoms in Eqs. © and (jU, since 
it appears in S only quadratic form: 

S = Sc -\- S f -\- Sj 

with 

Sc = T.{ll [^c(r, m,ir, t + A)- /i^,(r, t)vl/,(r, t)] - r J] ^,{r, t)^,{r', t)} 

t r (j.,r') 

^/ = E { E [^/('^^ ^)^/(^' t + A)- /i^^(r, t)mf{r, t)]]. 

t r 

The interaction between the two types of atoms is given by 

Si = U'J2^f{r,t)^f{r,t)^,{r,t)^,{r,t). 

r,t 

The integration over the Grassmann field in Z gives a space-diagonal determinant 

f e-^f-^'Y[d^f{r,t)d^f{r,t) = det{-dt + ft - t/'^'c^c), (12) 
where dt is the time-shift operator 

'*(r,t + A) A<t<(3 
.-^(r,A) t = P 

The second equation is a consequence of the antiperiodic boundary condition of the Grass- 
mann field. 

A. 2. Expansion with Ising Spins 

The partition function is now a functional integral of the c-Grassmann field 

Z = f e-''^deti-dt + ft-U'^,^c)V[^c]= f e-'^l[[l + l[m-U'^,ir,t)^,{r,t))]v[^,]. 
J J J. ^ 
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The product can be expanded in terms of Ising spins {S{r) = ±1} jl5[ as 



1 + S(r) 



l[[l + l[{fi-U'^,ir,t)^,{r,t))] = nn[/^-f^'^^(^'^)^c(r,t)] — 

r t {S{r)=±l} r t 

This reads with I = as 

{S(r)=±l} r 

Now the partition function Z can be expressed by a summation over configurations of Ising 
spins as 

Z= Z{{S{r)}) 

{S(r)=±l} 

with 

Z({5(r)}) = / e-^^ n /i^'^('-)e-(^'/M)iME.*=('-.*)*=(-.*)p[^^]. (13) 
J J. 

A. 3. Integration of the Light Atoms 

The c-Grassmann field appears only in a quadratic form in the partition function: 

After performing the \l/c^integration we obtain 

=/i^'S.iMdet(-9t + /i + £- — I). 

Following the same procedure for the Green's function, we obtain for G in Eq. (j3]) the matrix 

+ /i + t - (t/7/x)I)-i;t/^'S'-^(^)det(-9i + /i + £ - f I) 



G 
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